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ABSTRACT 

Neutrinos give a novel probe to explore deep interior of astrophysical objects, which otherwise is not 
accessible with optical observations; among notable examples are solar and supernova neutrinos. We 
show that there is a new class of strong neutrino emission from helium burning, 1 N + a — > 18 F + 7 
followed by beta decay 18 F — > 18 O + e + + u e , that gives a maximum neutrino luminosity of 10 8 times 
the solar bolometric luminosity at the helium-core flash of a 1 M star, whereas the flash is not 
observable by optical means. This means that the neutrino flux, of average energy of 0.382 MeV, will 
be 10% the solar CNO neutrino flux on Earth if the star is located at lOpc. 
Subject headings: stars: evolution — neutrinos 



1. PRELIMINARIES 

Low mass stars (M < 2M©) are known to develop a 
degenerate helium core during their e yolution along the 
Red Giant Branch (see for instance ISweigart fc GrossI 
dHI). The core is fed by the surrounding hydrogen- 
burning shell that slowly moves outwards in mass as 
hydrogen burning proceeds through the CNO bi-cycle. 
The helium mass fraction in the hydrogen-depleted core 
is Y = 1 — Z, where Z is the total heavy element mass 
fraction (note that Z does not change during hydrogen 
burning). As a result of the CNO burning, Z is com- 
posed by approximately 70% of 14 N. When the mass 
of the degenerate core reaches about 0.47M©, the tem- 
perature rises up to 10 8 K and helium ignites under de- 
generate conditions. At the moment of helium ignition, 
the total number of 14 N nuclei in the core is n( 14 N) w 
5.3 x 10 53 (M c /0.47) (Z/0.019), where M c is the core mass 
in solar mass units. At this temperature destruction of 
14 N during helium burning (independently of the stellar 
mass) occurs through the reaction 14 N + a — > 18 F + 7 
followed by beta decay 18 F -> 18 + e+ + v e . The 
maximum energy of these neutrinos is 0.633 MeV, while 
the average energy 0.382 MeV. This implies that were 
all the 14 N burned as described above, a total energy 
e v « 3.2 x 10 47 (M c /0.47) (Z/0.019) erg would be car- 
ried away from the star by neutrinos. The importance 
of this neutrino emission contrasts with the folklore that 
neutrino emission from nuclear reactions is not impor- 
tant during the helium burning phase. It is the purpose 
of this work to show, using stellar evolutionary calcula- 
tions, that most of this energy is actually released in the 
very short time-scale that characterizes the helium-core 
flash and that the neutrino flux reaches values as high as 
6 x 10 47 s" 1 . 

2. THE HELIUM-CORE FLASH 

The detailed structure of the degenerate cores of low 
mass RGB stars at the moment of h elium ignition has 
been widely discussed in the literature iSweigart fc GrossI 
119781: ISweigart et alJl9"89tlSalaris et al.l2002|) . Although 



uncertainties in the modeling exist, a very good under- 
standing of RGB stars has been achieved . This is re- 
flected, for instance, in the work bv lRaffeitl l)1990t) where 
it has been shown that observations of globular cluster 
stars and RR Lyrae stars indicate that the core mass 
at helium-flash ignition inferred from observations and 
those predicted by stellar evolution calculations agree at 
about a 5% level. This has been used to constrain differ- 
ent properties of neu trinos and also on some mor e exotic 
particles like axions i|Raffelt & Weissl|l992l Il995|) . 

In this work, however, we are particularly interested 
in nuclear neutrino production during the helium-core 
flash, which is a less well understood evolutionary phase. 
In particular, the central question that remains to be an- 
swered in a satisfactorily way is whether the helium-core 
flash is a dynamic event and whether dynamical episodes 
are of any consequence to the flash itself and the later 
evoluti on. However, hydrodynamic 2-D and 3-D calcu- 
lations l)Deupredl 19961: lAchatzlll995t and Achatz, Miiller 
& Weiss, in preparation) show that hydrostatic stellar 
evolution calculations appear to follow the global prop- 
erties of the helium-core flash reasonably well, at least in 
its global properties. With the confidence gained from 
those results, we use f or the present work the LPCODE 
stellar evolution code ( Alth aus et al]l2003lh Here we re- 
view some important aspects of particular relevance for 
the present work. The nuclear network has been up- 
date and now includes 33 isotopes from H to 32 P and 
99 r eactions (charged p articles reaction rates are taken 
from Angulo et al. 1999). Mixing in convective regions is 
treated as a diffusive process and equations arc fully cou- 
pled to nuclear burning. Conducti ve opacities are taken 
from iHubbard fc Lampardl l)1969|) . Neutrino emission 
rates are from Itoh and collaborators, exception made of 
the p lasma neutr i no em ission rate that has been adopted 
fromlHaft et all <I1994[) . More details can be found in 
lAlthaus et al.l l)2003|) . 

In what follows, we present the results for the evolu- 
tion of a 1 M Q stellar model of initial solar composition 
(Y,Z) = (0.275,0.019), our reference calculation. 
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Let us first briefly discuss the evolution of the model up 
to the helium-core flash. Calculations were started from 
a fully convective pre-main sequence model. The main 
sequence phase lasts about 9.5 Gyr and when hydrogen is 
depleted in the core the star evolves to the RGB phase, 
where the degenerate helium-core grows as burning in 
the hydrogen shell proceeds. As the mass of the helium 
core increases, it contracts and releases gravitational en- 
ergy that heats it up. Additionally, the temperature of 
the burning shell is mainly determined by the local value 
of the gravitational potential and thus increases steadily, 
also contributing to increase the core temperature. The 
dominating cooling mechanism in the core is the plasma 
neutrino emission, that also increases along with the core 
size. The net result is that the core temperature increases 
steadily along the RGB evolution, but neutrino emission 
is enhanced closer to the center in the almost isothermal 
core, eventually leading to a temperature inversion. In 
our evolutionary sequence, this happens when the degen- 
erate core has 0.27M Q and the age and luminosity of the 
model are 12.62 Gyr and log(L/L0)=1.8 respectively. 
The model continues to evolve along the RGB with its 
luminosity and core temperature increasing continuously. 
The core degeneracy also continues to increase, decou- 
pling the pressure from the temperature. This continues 
until the temperature in the core is finally enough to 
ignite the helium and the evolution along the RGB is 
halted. This defines the tip of the RGB. In our calcula- 
tions, the luminosity at the tip of the RGB for this model 
is log(-Ltip/Lo)=3.446 and the degenerate helium-core 
has 0.476 M . The age of the model is 12.664 Gyr. The 
temperature inversion is located at M max = 0.23M Q ; this 
sets the position where helium ignites. At the moment 
when the helium-burning luminosity Lhc exceeds for the 
first time that of the hydrogen-burning, the temperature 
T m ax and density /9 max at M max are 1.1 x 10 8 K and 
3.2 x 10 5 g cm" 3 respectively; the degeneracy parameter 
ijKippenhahn & WeieertJl 19901 Sect. 15.4) is ip = 8.5. 

Given the conditions at the moment of helium ignition, 
the core is supported by electronic degeneracy pressure 
and, as mentioned before, temperature is decoupled from 
pressure. The helium flash proceeds through the triple 
alpha (3a) reaction, that has a Maxwellian averaged 
rate (uv) 3a very strongly dependent on temperature: 

N l (<™>3a a r40 ""' 20 for T = X > -> 2 - 5 x 108 K - Tllis im- 
plies that a small increase in temperature results in an 

enormous increase in the burning rate, and because the 
structure of the star cannot respond to an increase in T 
by expanding and cooling, a thermonuclear runaway, the 
helium-core flash, quickly follows helium ignition. Time 
evolution during the helium-core flash of some important 
quantities is shown in Fig. ^ The top panel presents 
the time evolution of helium burning luminosity Lhc- In 
just a few years, L^ c grows from being a negligible en- 
ergy source in the star to overcome the luminosity of the 
star (Xboi) by almost seven orders of magnitude. It is 
important to note that when Lhc reaches about 10 2 L Q , 
energy cannot be transported by radiation alone, and 
adiabatic convection sets in in the layers above M max . 
At its maximum extension, the convective shell ranges 
from Mmax up to 0.475M Q . As shown in the next sec- 
tion, the extension of this convective shell determines to 
a large extent the total amount of N 14 that is actually 



burnt during the flash. Fig. also shows that the helium 
flash does not show up in Lboi, which actually decreases 
as the hydrogen burning luminosity (£h) drops abruptly 
when the layers surrounding the core expand and cool. 
To illustrate the conditions under which the flash starts 
and evolves, the middle and lower panels in Fig. ^ show 
the temperature and density at A/ max . The lower panel, 
in particular, makes evident that the energy release dur- 
ing the flash effectively removes the degeneracy in the 
burning region and is used to lift up the stellar core from 
its potential well. 

After the main episode of the helium-core flash, de- 
scribed above, evolutionary calculations show that a se- 
ries of sub-flashes follow, each of them occuring suces- 
sively closer to the center until degeneracy is completely 
removed from the core and s teady helium-burning is es- 
tablished |Serenell^^3^Iil l|2005f) . Detailed evolution 
through this phase is rather uncertain, in particular be- 
cause it is not well understood how heat f lows from the 
point of maxi mum temperature inwards (Achatz 1995; 
lDeupred ll996) . However, we do not expect that this af- 
fects our main results, because the bulk of the nuclear 
neutrinos and in particular the sharp peak in the neu- 
trino luminosity occurs during the main flash episode. 
We show this in the next section. 

3. NEUTRINOS FROM THE NITROGEN FLASH 

We now focus on the production of nuclear neutrinos 
during the helium-core flash. As mentioned before, 14 N 
is, by far, the second main consituent of the helium core 
at the time of helium flash ignition. The Maxwellian av- 
eraged rate of the 14 N + a — > 18 F + 7 reaction, Na (cv) 14 
(N^4 is the Avogadro number), is between 3 and 5 or- 
ders of magnitude larger than that of the 3a reaction 
in the temperature rang e characteristic of the helium 
flash l|Angulo et al.lll999j) and its dependence on tem- 
perature is extremely high (Na(o'v) 14 oc T 45 '-"' 21 for 
T = 1, ...,2.5 x 10 8 K). The linear dependence on p of 
the 14 N + a reaction delays its activation with respect to 
helium ignition through 3a, but the extreme temperature 
dependence leads to a very steep increase in the rate as 
temperature rises 1 . The 14 N half- life for a-captures de- 
creases from 10 15 s when helium ignites at T = 10 8 K 
to 1.7 x 10 4 s at T = 2.25 x 10 8 K, the maximum tem- 
perature reached in the flash. We call th is the "Nitro- 
gen Flash" , a term coined by llbenl l)1967l) in a pioneer- 
ing work on helium-flash calculations. At that time, the 
14 N + a reaction was thought to halt the evolution along 
the Red Giant Branch and actually lead to the ignition 
of the helium-flash. It seems clear now that, as described 
above, this is not the case and 14 N burning lags slightly 
behind the development of the flash, but we believe the 
name is still appropriate. 

The complete chain for 14 N burning is 
14 N + a ^ 18 F + 7 followed by beta decay 
i8p q _|_ e + _|_ v ^ where the neutrinos are pro- 
duced with a maximum energy of 0.633 MeV. The beta 

1 The helium flash is, at all times, driven by the 3« reaction, 
while 14 N burning supplies at most 10% of the nuclear energy re- 
lease and for a very short time. This is so due to a combination 
of factors: when the flash starts p > 10 5 g cm -3 and 3a domi- 
nates because its rate depends as p 2 . At later times, when the flash 
develops and starts producing n expansion of the burning region, 
density drops but 14 N is rapidly depleted. 
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decay half-life of 18 F is only 6.6 x 10 3 s, shorter than 
the half-life of 14 N for a-captures. Other destruction 
channels of 18 F are at most two orders of magnitude 
slower than the beta decay and do not play any sig- 
nificant role in the 18 F nucleosynthesis (in particular 
18 F(n, a) 15 N dominates among these channels when a 
small mass fraction of neutrons of about 10 -19 builds 
up transiently at the base of the burning shell). This 
has the consequence that the neutrino production 
rate basically follows that of 14 N destruction. In the 
upper panel of Fig. [2] we show the time evolution 
of the nuclear neutrino luminosity produced by 18 F 
decay. Note the very short time-scale that dominates 
its evolution, more than 4 orders of magnitude increase 
in only 0.2 years. The maximum luminosity is about 
10 8 L Q = 3.8 x 10 41 erg s _1 and the resulting number 
flux of neutrinos at the peak is about 6 x 10 47 s _1 . We 
have shown five epochs of the internal distribution of the 
14 N mass fraction in the lower panel of Fig. in order to 
illustrate how 14 N depletion and, accordingly, neutrino 
production proceeds. 14 N in the shell of M = 0.23 
to 0.475 M Q is destroyed nearly entirely within 1 yr. 
Indeed, 2.65 x 10~ 3 M Q of 14 N, which is 50% of 14 N 
produced in the previous CNO burning is depleted 
in this short time-scale. The nuclear neutrino energy 
emitted in nitrogen flash amounts to 0.05% the value 
emitted in the entire life of the star. 18 produced in 
the beta decay is eventually all consumed to produce 
22 Ne by further a-captures, which has a reaction rate 
comparable to the 14 N+a reaction, and therefore it 
entirely disappears in the core of white dwarfs. The only 
in-principle visible event that signatures helium flash 
are the burst neutrinos. 

The peak luminosity of nuclear neutrinos during the 
helium-core flash (LP cak ) is directly linked to the strength 
of the helium-flash and thus uncertainties in the helium- 
flash calculations will also affect the value of (XP cak ). A 
number of factors contribute to the uncertainty in the 
helium flash calculations. On one hand, uncertainties 
in the input physics is amply dominated by uncertain- 
ties of a factor of 2 in the electron conduc tive opacities 
(Kmn i i) in deg enerate cores of RGB stars l|Salaris et al.l 
2002r lCatela ll l2T)05|) . 

Other possible sources of uncertainty like thermal neu- 
trino emission, equation of state or radiative opacities are 
much better constrained and have a much s maller impact 
on th e uncertainties of RGB stellar models l|Salaris et alJ 
2002); and consequently in the nuclear neutrino emis- 
sion. A second source of uncertainty, that cannot cur- 
rently be quatified, is that introduced by the lack of 
complete hydrodynamic calculations of the helium-core 
flash. Although these calculations seems to support hy- 
dr odynamic stel lar evolutionary calculations (in particu- 
lar ^^h^£^^2^ shows that no mixing beyond the canon- 
ical convective boundaries as given by the Schwarzschild 
criterion occur and also that the temperature gradient in 
the convective shell is very close to the adiabatic one), 
uncertainties in how heat flows below the temperature 
inversion are present and can have some influence on the 
temporal evolution of the temperature at M max . Finally, 
the peak luminosity of nuclear neutrinos depends on the 
mass and composition of the star. We have addressed 
the first and last points by performing a series of ad- 
ditional evolutionary calculations. We do not attempt 



to cover the parameter space, but just to have a rough 
estimation of the level of the uncertainties that can be 
expected. The results are summarized in Tabled The 
behavior of Lf, cak at peak on mass and helium content 
reflects that of L^ e at peak. The behavior on the metal- 
licity follows approximately the linear relation expected 
from the simple arguments presented in Sect. ^ The 
dependence of Lf, cak on K con d seems to be linear, and a 
factor of 2 uncertainty in K con( j translates into a factor 
of 2 unceratainty in £P eak at peak. A last source of un- 
certainty, affecting only the neutrino production, is the 
N 14 + a reaction rate. Its unc ertainty direct l y tran slates 
into an uncertainty in iP cak . iGorres et alJ l)2000D have 
determined this reaction rate t o be 2 to 3 times smaller 
than that of lAngulo et al.ll999l (the value adopted in this 
work), with an uncertainty of about 20% in the temper- 
ature range of interest. 

4. CONCLUDING REMARKS 

We have shown that during the main episode of the 
helium-core flash a burst of neutrinos is produced due 
to the extremely rapid burning of N 14 by a captures 
in the so-called Nitrogen Flash. Our reference calcula- 
tion, a 1 M stellar model of solar composition (Y, Z) = 
(0.275, 0.019) gives a nuclear neutrino peak luminosity 
of about 1O 8 L . Neutrinos from the Nitrogen Flash 
carry precious information about the degenerate condi- 
tions under which the helium-core flash is believed to 
occur. They may be the only direct manifestation we 
can ever detect of such an event, that does not manifest 
itself in the stellar surface. The neutrino flux expected at 
the flash is 5.3 x 10 7 (o?/10pc) _2 cm _2 s _1 at maximum 
(for a star located d pc away), which lasts for 3 days. 
This is about 10% of the CNO neutrino flux from the 
Sun. 

We are not optimistic about the detection of these 
burst neutrinos, not only because this is a rare event, 
but also because solar neutrinos will give a serious back- 
ground. A very large volume detector with the capability 
of energy resolution with the threshold set just above the 
maximum energy of p-p neutrinos is needed. To give an 
idea, in a 1000 ton hydrocarbon scintillator, we expect 
for a helium flash at 10 pc 10 events for 3 days, which is 
to be compared to 80 events from solar CNO neutrinos 
for the neutrino energy window of 0.42 to 0.63 Mev (here 
neutrino oscillation is ignored). In practice, we also ex- 
pect a background from solar Be neutrinos for electron 
scattering. 

We want to acknowledge James Gunn and Bohdan 
Paczyhski for valuable discussions and comments, John 
Bahcall and Achim Weiss for a careful reading of the 
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Fig. 1. — Temporal evolution of stellar quantities during the 
development of the helium-core flash in our 1 MQStellar model. 
Top panel shows the evolution of helium (Lho) and hydrogen (Lh) 
burning luminosities, and stellar luminosity (Z/bol) m solar units 
(Lq= 3.84 X 10 33 erg s — 1 ). The stellar luminosity has been shifted 
up by 0.2 dex to make the otherwise two degenerate lines (Ljj 
and Lbol) visible. Middle and lower panels show, respectively, the 
temperature and density at the point of maximum temperature in 
the models (M max ), located off-center, at about 0.23 Mq. t = 
corresponds to the moment of maximum Lho- 



TABLE 1 

Important quantities characterizing the helium-core flash 
OF OUR evolutionary sequences. 



Model 



M e logL Ti p logL Hc log-Ly 



Reference 


0.476 


3.446 


9.84 


7.95 


0.5 X K con d 


0.484 


3.494 


10.18 


8.29 


2 x k coi1c j 


0.468 


3.403 


9.58 


7.65 


Y= 0.2505 


0.480 


3.440 


9.97 


8.12 


Y= 0.3005 


0.472 


3.430 


9.71 


7.78 


Z= 0.029 


0.474 


3.442 


9.78 


8.07 


Z= 0.009 


0.479 


3.413 


9.90 


7.70 


Z= 0.0019 


0.482 


3.347 


10.02 


7.14 


M= 0.80 


0.478 


3.432 


9.88 


8.01 


M= 1.20 


0.476 


3.434 


9.83 


7.94 


M= 1.40 


0.476 


3.434 


9.86 


7.95 



Note. — The table summarizes the most relevant quantities 
of our evolutionary sequences. First column denotes the model 
by giving the parameter that has been modified with respect to 
the reference model (first row). Second and third columns give 
the helium-core size and luminosity at the RGB-tip. Fourth and 
fifth columns give the peak helium-burning and nuclear neutrino 
luminosities during the flash. In all cases luminosities are in solar 
units. 
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Fig. 2. — Top panel. The evolution of the nuclear neutrinos 
luminosity (L v ) very close to the flash peak (t = 0). Note that 
nuclear neutrinos during the evolutionary phase shown here are 
solely created in the beta decay 18 F — > +e+ + v e that follows the 
o-capturc 14 N + a -> 18 F + 7. Vertical lines, labeled A, B, C, D 
and E denote five chosen epochs for which the internal abundance 
profile of 14 N in mass fraction is shown in the lower panel. Initially 
(A) 14 N is barely depleted. The temperature dependence of the 
reaction 14 N(a, 7) 18 F is even larger than that of the 3a reaction 
and when 14 N(a, 7) 18 F becomes fully active it produces a burst 
of neutrinos that lasts for less than 4 days (B-C) during which 
the neutrino luminosity of the star reaches a maximum value of 
about 1O 8 L0 (this is about a 5% percent of the total luminosity 
of nuclear neutrinos of the Milky Way). 50% of the available 14 N 
in the burning region is burnt in this time interval. Strong 14 N 
burning continues (D) but removal of degeneracy and subsequent 
expansion and cooling of the burning layers quickly follow and set 
a much slower pace for 1 N depiction (E). The step-like shape 
of the profiles denotes the presence of convection, that keeps the 
shell with an almost homogeneous composition even when nuclear 
burning is actually occuring only in the 0.23 — O.27M0 region. 
In passing by, we note that the shape of the 14 N profile in the 
inner 0.2MQreflects the increasing temperature of hydrogen-shell 
burning with increasing core size during previous evolutionary 
stages. 



